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The wedge geometry closed by a circular-cylindrical arc is a nontrivial gen- 
eralization of the cylinder, which may have various applications. If the radial 
boundaries are not perfect conductors, the angular eigenvalues are only im- 
plicitly determined. When the speed of light is the same on both sides of the 
wedge, the Casimir energy is finite, unlike the case of a perfect conductor, 
where there is a divergence associated with the corners where the radial planes 
meet the circular arc. We advance the study of this system by reporting results 
on the temperature dependence for the conducting situation. We also discuss 
the appropriate choice of the electromagnetic energy- momentum tensor. 



1. Introduction 

Casimir theory for the wedge geometry continues to attract interest. The 
reasons for this are many-faceted - probably the most important one be- 
ing that the material boundaries are plane, thus avoiding some of the for- 
mal divergences that so often plague calculations in the presence of curved 
boundaries. The wedge geometry moreover implies a formalism closely re- 
lated to that of cylindrical geometry, and actually also to that of cosmic 
string theory. Finally, the wedge geometry is a convenient testing ground 
for experimental tests of Casimir-Polder forces. 

The Casimir energy and stress in a wedge geometry was approached 
already in the 1970s. 1,2 Various embodiments of the wedge with perfectly 
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conducting walls were treated by Brevik and co-workers 3-5 and others. 6,7 
More recently a wedge intercut by a cylindrical shell was considered by 
Ncstcrcnko and co-workers. 8 ' 9 Local Casimir stresses were considered by 
Saharian and co-workers. 10-12 The interaction of an atom with a wedge was 
studied experimentally by Sukenik et al. 13 The theory of that interaction 
was worked out by Barton 14 and others. 15-18 The semitransparent wedge 
has very recently been considered by Milton, Wagner, and Kirsten. 19 The 
closely related case of circular symmetry has been treated in several papers, 
dealing with a perfectly conducting circular boundary, 20-22 as well as the 
case of a dielectric circular boundary. 24-27 

The typical wedge geometry is sketched in figure la. The planes are 
situated at 8 = and 9 = a. We shall assume in the present paper that the 
interior of the wedge is filled with an isotropic medium of spatially constant 
and nondispersive refractive index n = ^/sjl. 



(a) (b) 




Fig. 1. Wedge geometries, (a) The perfectly conducting wedge geometry, (b) The ge- 
ometry of a wedge intercut by a perfectly conducting cylindrical arc. (c) Wedge with 
magnetodielectric arc. (d) Diaphanous wedge in a perfectly conducting cylindrical shell. 

In the simplest version of the wedge model, the planes are taken to be 
perfectly conducting. Various modifications of this simple wedge model can 
be envisaged. In sections 3 and 4 below we consider two generalizations 
of the simple wedge geometry, called Wedge I and II, in which the interior 
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region is closed by a circular boundary thus implying an eigenvalue problem 
for the photon frequencies. The Wedge II model treated in section 4, in 
particular, removes the strict perfect boundary condition of the radial walls. 
The material of these two sections is based on two recent papers. 28,29 New 
developments are a closer examination of the behavior at finite temperature. 
As an introductory step, we delineate in the next section the essentials 
of classic Casimir theory for the perfectly conducting wedge. In section 6 
we relate the macroscopic electromagnetic theory for the dielectric wedge 
region to the more general question about which electromagnetic energy- 
momentum tensor is to be preferred in media. This one-hundred-year-old 
question has actually attracted considerable interest recently. 

2. Extracts from the Classic Theory for the Perfectly 
Conducting Wedge 

The geometry is shown in Fig. la and was considered in a number of pub- 
lications. 1-5 The governing equation for the Fourier transform of Green's 
function T(x, x') is 

V x V x T(r, r', w) - e^ 2 r(r, r', lu) = -fiuj 2 15(r - r'). (1) 

After solving this equation in terms of the scalar Green's functions F m (r, r') 
and G m (r,r') 5 we can calculate the effective field products for the electric 
fields as i(E i (r)Ek(r') u) — r^r, r', u). The corresponding products for the 
magnetic fields follow from Maxwell's equations. The points r and r' are 
assumed to be close but not coincident. 

The effective products can now be inserted in the electromagnetic 
energy- momentum tensor whose spatial part in classical notation is, 
in Heaviside-Lorentz units, 

S ik = -EiD k - H t B k + U lk (V D + H B). (2) 

A lengthy calculation leads to the expression (Minkowski metric assumed) 

( VM> = (J + ll) - l) diag(l, -3, 1, 1), (3) 

where we have subtracted off the term {S"^ w (r)) corresponding to a plane 
sheet. The components are ordered as {S^} = (S rr , See, S zz ,—w), where 
w is the energy density. 

The expression (3) refers to zero temperature. It is worth noticing that 
in the limit a — > 0, r — > oo such that ar becomes the separation between 
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parallel plates, the expression agrees with Barton. 1 Moreover, the expres- 
sion agrees with that of a cosmic string if the string's deficit angle 8nGM 
is identified with 2ir - 2a. 28 

3. Wedge I: Perfectly Conducting Walls; Circular 
Boundary at r — a 

The geometry is shown in Fig. lb and lc. The walls are perfectly conducting 
as before, while we assume now that there is as boundary a circular arc 
with radius a. The wedge thus has an interior region r < a (region 1) 
where the material parameters are e\ and /ii, and an exterior region r > a 
(region 2) with analogous parameters and fi2- The materials are assumed 
to be nondispersive. We assume the diaphanous (isorcfractive) condition 
eifii = £2M2 = n 2 - The transverse wave numbers k± in the two regions 
are accordingly the same, k\ = n 2 uj 2 — k 2 . It is convenient to introduce 
the symbol p = n/a, and also A„(:r) = (I^(x)K u (x))' ,where I v and K v are 
modified Bessel functions. 

3.1. The boundary r — a being perfectly conducting 

This is the simplest case. Detailed expansions of the electric and magnetic 
fields are given in Ref. 28. In region 1 there are two independent polar- 
izations, one TM polarization where the mode eigenvalues are determined 
by J mp {k±a) — with m = 1,2,3, and one TE polarization where the 
eigenvalues are determined by J' mp (k^_a) = 0. One azimuthally symmetric 
TE mode exists, corresponding to m = 0, but there is no such TM mode. 
In region 2 the TM polarization yields Hm P {k±a) = 0, m = 1,2,3,.., 

whereas the TE polarization yields Hmp (k±a) = 0, m = 0, 1, 2, Sum- 
ming over all modes and making use of the argument principle, we arrive 
at the following expression for the total zero-point energy 



the prime meaning that the mode m = is counted with half weight. 
This is the boundary-induced contribution to the zero-point energy. If the 
boundary r = a were removed and cither the interior or the exterior medium 
were to fill the whole region, we would get £ = 0. Moreover, we have omitted 
a zero-mode divergence caused by the sharp corners where the arc meets the 
wedge. If p — 1, Eq. (4) is one-half that for a conducting circular cylinder. 
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3.2. Dielectric boundary at r — a 

The most important change compared with the previous subsection is that 
regions 1 and 2 become coupled via electromagnetic boundary conditions. 
Assuming m = n 2 we hnd the following simple expression for the zero-point 
energy 

*=X=^E'/ <fc*Ml-fVAU (5) 



Airna 2 / — ' 

m=0 



with 



£=^^. (6) 

£2 + ei K ' 

The conducting case is obtained by setting £ = 1. The case n\ ^ n 2 is more 
complicated, but for weak-coupling, \e\ — e 2 \ «l,a self-energy can still be 
extracted 28 by generalizing the work done on dielectric cylinders. 24,25 

The expressions (4) and (5) are still not in general finite. A finite self- 
energy can be extracted from this formula by a method of zeta function reg- 
ularisation, 23 generalizing the standard formal result for a circular cylinder 
(p = l). 20 Further details and numerical results are reported in Ref. 28. 



4. Wedge II: Diaphanous (Isorefractive) Wedge in 
cylindrical shell. 

Consider the geometry of Fig. Id wherein a diaphanous magnctodiclcctric 
wedge (rii — n 2 ) inside a perfectly conducting cylindrical shell of radius a 
is considered. More details were published in Ref. 29. The sum over orders 
of the Bessel function partial waves is now not simply equidistant values 
v = mp as before but the zeros of the dispersion function for v = in: 

D(ir); a) — sinh 2 r/ir — £ 2 sinh 2 t](t: — a). (7) 

The reflection coefficient £ was defined in Eq. (6). In the absence of any 
wedge this becomes Do(irj) = sinh 2 nir. The energy of the diaphanous wedge 
enclosed by a perfectly conducting cylindrical shell is thus found as 

* - a* £ * £ d «£ 11,11 - ^ (8) 

In the non-dispersive case where £ is independent of £ this may be 
simplified by means of partial integration w.r.t. £, introduction of polar 
coordinates, and use of symmetry properties, to 
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Since Re{arAf„} < 1 we may use for numerical purposes 



„ „ x 2 ImA? 
Imln(l - x 2 \j ' ) = - arctan 



l-^ReA? n - 
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Fig. 2. The energy per length of the wedge, divided by a reference energy (8irna 2 ) 1 , 
as a function of opening angle a for different values of the reflection coefficient f . 



A numerical evaluation of this energy expression was performed in 
Ref. 29 and the details of the procedure will not be iterated here. The 
energy is plotted as a function of the opening angle a for different values 
of the reflection coefficient |£| in Fig. 2. 

5. Considerations of finite temperature 

We present for the first time some considerations on the Casimir energy of 
a closed perfectly conducting wedge such as in figure lb when T > 0. (For 
earlier work on the cylinder at high temperatures see Ref. 30.) By letting 
p = 1 in the following, these considerations automatically apply to the 
case of the perfectly conducting cylinder. We consider the finite part of the 
energy, given for zero T in equation (4). As is customary, finite temperature 
implies a compactification of the imaginary time axis so that the integral 
over imaginary frequencies i£ becomes a sum over Matsubara frequencies 



dC/(C)->2^T^'/(0). 



(10) 
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A little calculation shows that the expression (4) may then be written 
£ = J2m'=o £ "i where 



^ = ^E' r , I™-* 2 HI ~ ^(')1 = ^ E W (11) 

The nondimensional temperature is r = 2irnaT. As previously encountered 
at zero temperature the simple expression for £ is divergent and the chal- 
lenge is to regularize it. We follow a scheme closely analogous to that of 
Ref. 23 (c.f. appendix A of Ref. 28) using zeta functions by subtracting and 
adding the asymptotic behaviour of the integrand, 

x 4 

ln[l - x 2 X 2 mp (x)\ ~ - 4(m2p2 + ^3 , m,x^oo. (12) 

This asymptotic behaviour is responsible for divergences both for j — > oo 
and m — > oo. The work on this problem is still in progress and details of 
the calculations and further discussion will be published elsewhere. 31 
Adding and subtracting the asymptotic behaviour we write 



rp OO 

£ =—^' 



ira 



3=0 L 



; 1 f°° dxx 5 
4 J jT 



a/ x 2 — j 2 T 2 (m 2 p 2 + x 2 ) 3 



(13) 



using compact notation m defined as m(m) = m for m > 1 and m(0) = 1 /p. 
The symbol e implies that the leading asymptotic term (12) has been sub- 
tracted from the integrand, with m replacing m. The double sum resulting 
from the integral in (13) is formally divergent but may be regularised by 
use of the Chowla-Selberg formula. 32 After some calculation we obtain the 
following result 

,T fw. r(3-3rd T +r 2 d 2 )[l + 2JC(r)+S(r,p)} ^ 

m,j—0 
~i JT^l+J 2 T 2 

2 8 °° 

5(r,p) =-[7-ln(47rp/r)] + - V ct q {1)K q (2-kIt /p), (14c) 
P Pf~[ 

where 7 is Euler's constant, (7q(1) is the number of positive divisors of I, K 
is the modified Bessel function of the second kind, and d T = d/dr. Notably 
S obeys the symmetry relation S(r,p) — S(p, r). One may show that this 
result reduces to the zero temperature expression 28 when r — ► 0, 31 and we 
have verified that when p = 1 the coefficients of the two leading-order terms 



December 3, 2009 21:28 



WSPC - Proceedings Trim Size: 9in x 6in QFEXT09wedgeE 



8 

as t — ► oo , of order r and rlnr, equal twice those found previously for the 
cylindrical shell in Ref. 30 as they should. 

6. On the Electromagnetic Energy Momentum Tensor in 
Media 

The expression (2) for the spatial part of the energy-momentum tensor 
(equal to minus the Maxwell stress tensor according to M0ller 33 ), is 
common for the Minkowski and Abraham tensor alternatives as long as the 
medium is isotropic. The extraction of the "correct" form of the tensor has 
however been discussed for a long time. Thus it is to be noted that 

• the expression is different from that of Einstein and Laub (1908); 34 

• it is different from that of Peierls (1976); 35 

• and it is different from that of Raabe and Welsch (2005) ; 36 cf. also 
the comment 37 of Brevik and Ellingsen. A survey up to 1979 is given by 
Brevik. 38 

When combined with the Minkowski momentum density g M = D x B 
one obtains the Minkowski energy-momentum tensor S^f v whose covariant 
form can be written as 33 

Sffv = F^aH ua — -5^F a fjH a/ 3. (15) 

Here F 4k = iE k ,H 4k = iD k ,F lk = B t (cycl), H ik = Hi (cycl). This 
energy-momentum tensor has several attractive properties: it has zero di- 
vergence for a pure radiation field, the four components of momentum and 
energy forming a four-vector; it is a convenient expression for field quanti- 
zation (for instance, the quantization for a radiation field in a nondispersive 
medium can be found via a mapping technique leading one from vacuum 
to a medium 39 ), and it has the peculiar property of being a space-like ex- 
pression leading to negative photon energies in certain coordinate systems, 
strikingly found in connection with the Cherenkov effect in the emitter's 
rest frame. 

One might think that it should be relatively easy to test the Minkowski 
tensor by measuring electromagnetic forces in optics. Actually, this is not 
so easy as most experiments measure only the surface force density f = 
— 7}E 2 Ve, acting on surfaces. Let us give some examples: 

The classic experiment of Ashkin and Dziedzic 40 showed how a narrow 
light beam incident on a free liquid surface acts by giving rise to an outward 
pull. A related experiment is that of Zhang and Chang, 41 demonstrating 
the oscillations of a water droplet by a laser pulse. The optical stretcher 
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experiment of Guck et al. 42 is of the same kind, as is the series of two- 
fluid experiments of Dclvillc et al. 43 near the critical point. The important 
factor in these experiments is simply the surface force, not electromagnetic 
momentum. And the recent fiber experiment of She et al. 44 belongs in our 
opinion to the same category. 45 The experiment of Campbell et al. 4e is 
however different in nature, as it is one of the very few experiments being 
able to test the Minkowski momentum directly. 

Under stationary conditions, where the high-frequency forces average 
out when averaged over a period, the Minkowski tensor is able to describe all 
experiments that we are aware of. And this gives support to our expression 
(2) for the stress tensor. 

The work of KAM was supported by the US National Science Founda- 
tion under Grant No. PHY-0554926 and by the US Department of Energy 
under Grants Nos. DE-FG02-04ER41305 and DE-FG02-04ER-46140. 
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